Ordering kinetics in an fee A^B binary alloy model: Monte Carlo studies. 
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Using an atom-vacancy exchange algorithm, we investigate the kinetics of the order-disorder tran- 
sition in an fee A^B binary alloy model following a temperature quench from the disordered phase. 
We observe two clearly distinct ordering scenarios depending on whether the final temperature Tf 
falls above or below the ordering spinodal Tsp, which is deduced from simulations at equilibrium. 
For shallow quenches {Tf > Tsp) we identify an incubation time nnc which characterizes the onset 
of ordering through the formation of overcritical ordered nuclei. The algorithm we use together with 
experimental information on tracer diffusion in CusAu alloys allows us to estimate the physical time 
scale connected with Ti„c in that material. Deep quenches, Tf < Tsp, result in spinodal ordering. 
Coarsening processes at long times proceed substantially slower than predicted by the Lifshitz- 
Allen-Cahn t^^^ law. Structure factors related to the geometry of the two types of domain walls 
that appear in our model are found to be consistent with Porod's law in one and two dimensions. 

PACS numbers: 05.50.+q, 64.60.-1, 64.60.Cn 
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I. INTRODUCTION 

Several face-centered cubic binary metallic alloys, like 
CusAu, CuaPd, Mgsin, CoaPt, etc. exhibit long range 
order with a Ll2-structure below some specific ordering 
temperature Tq. In the Ll2-structure four equivalent 
phases exist, where one of the four simple cubic sublat- 
tices building the fcc-lattice are preferentially occupied 
by the minority atoms. In studying kinetic processes of 
phase ordering, the following general considerations must 
be taken into account: 

(i) The transition is of first order. Hence, for small 
enough supercoolings the disordered phase remains 
metastable. Relaxation at short times after a tem- 
perature quench is governed by the formation of 
ordered nuclei, which grow into the disordered ma- 
trix. On the other hand, the. concept of spinodal 
ordering has been advancedtl to characterize the 
phase ordering processes under large supercoolings. 

(a) The degeneracy of the ordered phase implies a four- 
component order parameter, ^ — (■0o, ■02, V's)- 
Here, ■j/'q with a = 1,2,3, are non-conserved struc- 
tural order parameter components coupled to a 
conserved density, ipo, which describes the concen- 
tration of the two atomic species. 

{Hi) The antiphase domain structure is anisotropic as a 
result of the existence of two types of antiphase do- 
main walls: low-energy (type-I) and high-energy 
(type-II) walls. 

Under these conditions a rich spectrum of kinetic pro- 
cesses is to be expected. In particular, there remain still 
open questions concerning scaling and universality in the 
late-stage growth kinetics. This has motivated several 
experimental investigations of ordering after a tempera- 
ture quench, mostly on CusAu, where direct information 



has been obtained from time-resolved X-ray diffraction 
and neutron scattering. On the other hand, only few 
theoretical or computer simulatiniLiitudies on such^mate- 
rials have been reported so far nua Frontera et alB have 
recently simulated the growth of Ll2-ordered domains 
both within an atom-atom exchange and the more realis- 
tic atom-vacancy exchange mechanisms, for quench tem- 
peratures below the expected spinodal temperature Tgp. 
A similar model with atom-vacanoy exchange {ABV- 
model) was applied by Kessler et al.a to surface-induced 
kinetic ordering processes near Cu3Au(001). 

In this paper we investigate the ordering kinetics in the 
bulk of such alloys at quench temperatures Tf both below 
and above Tgp. In contrast to earlier|-Monte Carlo stud- 
ies of nucleation in Ising-type modelsErlj we again employ 
the ABV~mode\ with effective chemical interactions be- 
tween nearest neighbors on the fcc-lattice. This allows 
us to compare our results directly with experiments on 
fcc-alloys, where vacancy-driven processes prevail. In 
particular, we analyze the nucleation regime, where we 
find a well defined incubation time that sensitively de- 
pends on the depth of the quench, Tq — Tj. Moreover, for 
Tf < Tsp, with Tsp deduced from static correlations, we 
observe domain patterns typical to spinodal ordering. In 
the long-time coarsening regime we extract anisotropic 
scaling functions. Our model is a limiting case where 
type-I domain walls have exactly zero formation-energy 
and therefore are extremely stable. Within our acces- 
sible time window we find coarsening exponents which 
are significantly smaller than the convenjtkinal exponent 
a = 1/2 for curvature-driven coarsening Eu 

After shortly explaining our model and simulation 
techniques in section II, we compare in section [II order- 



ing processes that occur for shallow (To > Tf > Tsp) and 
deep {Tf < Tsp) quenches. In section IV we define the 
incubation time Tine and investigate its dependence on 
the depth of the quench. Anisotropic scaling functions 
are discussed in section 0, while section VI contains a 
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short summary of our results. 

II. MODEL AND SIMULATION METHOD 

We consider a three-dimensional lattice oi LxLxL fcc- 
cells with cubic lattice constant a and periodic boundary 
conditions in all directions. Unless otherwise stated, we 
use L = 128. Each site, i, of the lattice is occupied either 
by an atom of type A, an atom of type _B, or a vacancy, 
with an obvious condition that the corresponding occupa- 
tion numbers fulfill cf + cf + — I- In accord with the 
stoichiometry of the A^B alloys there are exactly three 
times as many A atoms present as B atoms. The number 
of vacancies is chosen small enough so that they do not 
affect static properties of the system. 

In our simplified modela only nearest neighbor atom- 
atom interactions are taken into account. The corre- 
sponding Hamiltonian then reads: 



in terms of weighted averages of the occupation numbers 
computed of all sites belonging to the cell. (The weights 
are taken as the inverse of the numbers of neighboring 
fcc-cells that "share" the site in question, i. e. functions 
1/2 and 1/8 are taken for the "face" and "corner" sites, 
respectively) . 

As can be seen from Eq. (2), each of the non-conserved 
components of the order parameter describes a modula- 
tion of the i3-atom concentration along one of the cu- 
bic axes. For example, 7^ means that the system 
contains alternating B-enhanced and i3-depleted atomic 
layers along the x-axes. As a consequence of such layer- 
wise arrangement of B-atoms, the (100) peak shows up in 
an X-ray diffraction experiment in addition to the (200), 
(020), (002) and (111) peaks characteristic of the un- 
derlying fcc-lattice. Similarly, 7/^2 and ip^-type ordering 
leads to additional (010) and (001) peaks, respectively.EJ 
For a distance k from these superstructure peaks the scat- 
tering intensity is described by the structure factors. 



H 



[VAAcfcf + VBBcfcf + VABicfc 



B„B 
j 



cf)], 



^a(fc,t) = (|^'a(fc,t)|'). 



(3) 



(1) 



where the sum is restricted to nearest-neighbor pairs. 
We are interested here in the transition to the ordered 
Ll2-structure in which minority atoms [B) predomi- 
nantly occupy one of the simple-cubic sublattices of the 
original fcc-lattice. It is then natural to assume Vbb > 0, 
Vaa < and Vab < 0. The ground state of the 
LI2 phase is fourfold degenerate and corresponds to all 
B-atoms segregated to exactly one of the simple cubic 
sublattices. For a stoichiometric A^B-alloy the tran- 
sition occurs at a temperature Tq ~ 1.83|J|/A:b, with 
J = —jiVAA + Vbb — ^Vab) as discussed in Ref. ^ (cf. 
also Refs. p|,ppC|). The ordered phase is characterized 
by one conserved, scalar order parameter, -00, related to 
the composition and three non-conserved order param- 
eter components (V'l, V'2, V's)- These are defined by the 
following equations: 



ipo = mi + m2 + + m4, 

Ipl — mi - 771 2 - 77l3 + TO4, 
0^2 = mi - ?7l2 + 7713 - W4, 



(2) 



^3 



mi + m2 — — 



where 77ict are the differences of the mean A~ and 
B-occupation numbers, ma = (cf) — (cf } for i e a. 
The index a = 1 ... 4 enumerates the four equivalent 
simple-cubic sublattices of the fcc-structure. For a 
homogeneous A^B alloy, ipo — 2. In the disordered 
phase, 4'i — — 'ips = 0, while the four equiva- 
lent ordered phases are described as {ipi,ip2,''p3)/ip = 
(-1,1,1); (1,-1,1); (1,1^-1) and (-1,-1,-1), with a 
temperature-dependent (V' = 2 at T — 0). For a non- 
homogeneous system Eq. (2) gives the corresponding lo- 
cal order parameters for each (cubic) elementary fcc-cell 



which in a non-equilibrium state depend on time t. In 
Eq. (3), a{kx,ky,kz,t) denotes the Fourier transform 
of the order parameter ipaix, z,t), and a = 1,2,3. 
The widths of these intensity profiles are determined 
by the sizes of antiphase domains.EiJ As mentioned in 
the Introduction, this system displays two types of an- 
tiphase boundaries, low-energy type-I boundaries with 
no change in the arrangement ofoifiarest neighbors, and 
high-energy type-II boundaries.E3 Since the Hamilto- 
nian, Eq. (1), involves only nearest neighboipiiteractions, 
type-I boundaries in fact have zero energy.Ej 

We assume here that the time-evolution of our model 
occurs only through atom-vacancy exchange processes. 
In an elementary move we first randomly single out one 
of the vacancies and, second, we choose at random one of 
its nearest neighbors occupied by an atom. The standard 
Metropolis algorithm is then used to decide whether the 
atom-vacancy exchange for the chosen pair takes place 
or not. In one Monte Carlo step (MCS), the system com- 
pletes a series of as many such elementary moves as there 
are lattice sites in the system, i. e. 4L^, so that this time 
unit does not depend on jthe actual number of vacancies. 

In previous simulationsQ we have chosen the interaction 
parameters such that they were consistent with the ob- 
served Au-segregation at the CusAu (OOl)-surface. Fol- 
lowing Ref. 1^, we take Vbb — —Vaa — —Vab > and 
assume ~ 6.1 • 10~^|-as the mean density of vacancies 
present in the systemJl^ For vacancy concentrations of 
that order we have verified that both static and kinetic 
properties are independent of the precise value of . To 
simulate a sudden quench from high temperatures to a 
final temperature Tf < To we start with a random distri- 
bution of atoms at time t — and let the system evolve in 
time at Tf. The ensuing equilibration process is analyzed 
by calculating energies, structure factors, and other or- 
dering characteristics from averages over 10 independent 
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FIG. 1; Tracer diffusion coefficients of A (full squares) and 
B (full triangles) atoms in an A3B fcc-alloy. Physical units 
were assigned to the Monte Carlo data with the help of the 
experimentally known order-disorder transition temperature, 
To = 663-K, and temperature-dependent concentrations of 
vacanciestj in CusAu. The dasted line represents measured 
An tracer diffusion coefficients Ji3 The open circle marks the 
experimental point used to set the time scale. In the inset the 
corresponding "raw" Monte Carlo results are displayed. 



Monte Carlo runs. 

As the first application of our atom-vacancy exchange 
algorithm we have calculated tracer diffusion coefficients 
for A- and i?-atoms in the disordered phase close to 
Tq. Since in our model Vaa — Vab, the jump of an 
A-atom will leave the interaction energy with its envi- 
ronment almost unchanged and hence Da is found to be 
nearly temperature-independent. The i?_B-repulsion, on 
the other hand, gives rise to a temperature-dependent 
Db < Da, see the inset of Fig. 1. As these results refer 
to a fixed vacancy concentration, they cannot be directly 
related to experiments, where (T) normally is a strong 
function of T.Eil Nevertheless, the ratio Da/ Db ~ 2 at 
T/Tq = 1.7 very roughly agrees with the experimental 
valuCpDcuZ-Dyiji — 1-45 for CusAu at the same T/Tq- 
ratioJ13 thus supporting our choice of nearest-neighbor 
interaction parameters. Moreover, in an attempt to map 
the Monte Carlo time to the physical time scale, we can. 
exploit experimental knowledge of c^(r) for CusAu.tfl 
As <C 1, hopping rates of a tracer atom are simply 



proportional to 



This suggests to introduce a new 



Monte Carlo time unit 1 MCS* = (T) ■ 1 MCS, where 
1 MCS was defined above as one attempted vacancy ex- 
change per lattice site. Regarding the simulated mean- 
square displacements as a function of the new time scale 
(with units MCS*) allows us to extract diffusion con- 
stants Da and Db as shown in Fig. 1. Scales on the 
axes were obtained using the CusAu transition tempera- 
ture To — 663K, the lattice parameter a = 3.74A and an 
additional parameter, which is determined by equating 
the calculated diffusion coastant Db with the experimen- 
tal Dau at T/Tq = 1.27.li3 The last parameter converts 
1 MCS* directly into seconds. The favorable compari- 



son between calculated and measured diffusion constants 
suggests that the above procedure, based on the atom- 
vacancy exchange algorithm and the knowledge of the 
experimental c^(r)-function, provides a description of 
processes on the physical time scale. We shall come back 
to this issue in section III in the context of nucleation 
processes. 



III. THERMALLY ACTIVATED VERSUS 
CONTINUOUS ORDERING 

At r = To the model-system defined in the 
ing section undergoes a first order phase transitio: 
from the disordered phase (for T > Tq) to the Ll2-type 
ordered structure (for T < Tq). The correlation length 
of the disordered phase, ^, remains finite at coexistence, 
and approximately fulfills the mean-field type relation 
^^(rl oc 1/(T — Tsp), with Tsp < Tq. In our previous 



workd we found Tgp = (0.967 ±0.003)To as the value best 
fitting our Monte Carlo data for the correlation length 
above Tq. Hence, within this procedure, our model dis- 
plays a fairly well-defined temperature for the onset of 
spinodal ordering. A precise separation of the metastable 
(thermally activated ordering) from the unstable regime 
(continuous ordering) , howevet,pdpes not exist in systems 
with short-range interactions.BEZi 

Below we demonstrate the role of that tempera- 
ture deduced from equilibrium considerations, in non- 
equilibrium ordering processes. It turns out that for 
temperatures Tf > Tgp and T/ < Tsp two contrasting 
transition scenarios occur that lead to domain patterns 
typical of nucleation and spinodal ordering, respectively. 
We shall illustrate this by one example for each regime. 

Let us start our description with thermally activated 
growth of the ordered phase, i. e. with the case of a 
shallow quench, Tq > Tf > Tgp. We then expect that 
the ordered regions - nuclei of the new phase - are re- 
peatedly formed and destroyed by thermal fluctuations 
unless a nucleus exceeding a certain critical size is built. 
Once formed, such an overcritical nucleus will grow rel- 
atively fast against the surrounding disordered bulk un- 
til it meets another (growing) ordered domain. An xy- 
section of a configuration with overcritical nuclei formed 
in a system of 128 x 128 x 128 fcc-cells is shown in 
Figs. 2a and 2b. It emerged 7000 MCS after the quench 
to Tf = 0.972 To > Tsp. The central nucleus seen in 
the figure is in fact composed of three "twinned" crys- 
tallites, corresponding to i/ji ~ 2 (the large, nearly ho- 
mogeneous, bright spot in Fig. 2a) and tpi ~ —2 (the 
neighboring dark spot). As shown in Fig. 2b, the ho- 
mogeneous bright spot of Fig. 2a is in fact composed of 
two crystallites corresponding to -03 c± —2 and 7/13 ~ 2, 
respectively. A second nucleus characterized by i/'i — 2 
and -03 ~ —2 grows at the right edge of the system. The 
configuration displayed in Figure 2 is quite typical for a 
quench temperature slightly above Tsp. 

To follow the ordering process that takes place in the 
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a) Vi Tf = 0.972 To h) a) Vi Tf= 0.5465 Tg h) 




FIG. 2; Thermally activated ordering in a system of 128 x 
128 X 128 fcc-cells at T/ = 0.972 To > T,p. (a) ipi and 
(b) tp3 patterns containing overcritical nuclei at a time t = 
7 X 10"^ MCS after the quench from random initial condi- 
tions (grey scales indicate local values of t/ic (a = 1,3) be- 
tween —4 and 4). (c) Histograms of corresponding $ = 
+ \ip2 \ + \tp3\)/3 values averaged over blocs of 4 x 4 x 4 
fcc-cells calculated at a series of times after the quench. 



whole volume of the system we now introduce an addi- 
tional quantityEl $ = + |V'2| + |V'3|)/3 as a conve- 
nient indicator of the degree of local order: <i> ~ in the 
disordered phase and (f> < 2 for any of the four equiva- 
lent ordered states. Further, we divide the system into 
blocs of 4 X 4 X 4 fcc-cells, calculate $ independently for 
each bloc to obtain a distribution function P($,t) for a 
series of time intervals after the quench. For a quench 
to Tf — 0.972 To the obtained histograms are shown in 
Fig. 2c. In this figure we easily identify three stages of 
the ordering process. First, up to several thousand MCS, 
the whole system remains disordered. The distribution 
of $ broadens with time but stays concentrated close to 
$ = 0. Then, after some incubation time. Tine, see sec- 
tion 1^, the second maximum close to $ = 2 shows up. 
At that moment first overcritical nuclei emerge in the 
system and start to grow. The second stage, in which 
overcritical nuclei grow against the disordered bulk, is 
realized by histograms for 7 • 10^ < t < 10"* MCS in 
which both peaks are clearly visible. In the third stage. 



FIG. 3: Continuous ordering in a system of 128 x 128 x 128 fcc- 
cells at Tf = 0.5465 To < Tsp. (a) tpi and (b) Va patterns at 
a time t — 7 x 10'^ MCS after the quench from random initial 
conditions. Two types of domain walls are clearly shown, (c) 
Histograms of corresponding $ values calculated at a series 
of times after the quench. 



practically the whole system is already ordered and, as 
illustrated by curves for i = 2 • 10^ and 5 • 10"* MCS, the 
distribution of <i> changes very little with time. At that 
stage coarsening processes take place. 

In contrast to the above sequence of events, deep 
quenches (Tf < Tsp) result in continuous ordering, in 
which shortly after the quench the whole volume of the 
system decomposes into ordered domains. As a represen- 
tative example Figs. 3a and 3b display distributions of V'l 
and Tps, respectively, for an xy-section through the sys- 
tem at t = 7 • 10^ MCS after a quench to Tf = 0.5465 Tq. 
The two-tmes of domain walls, that may be formed in the 
system,aE3 show up in Figs. 3a and 3b. The low-energy 
interfaces appear there as straight lines parallel to the x 
or y-axis. Across each of these lines tpi and -03, or -02 and 
"03 simultaneously change sign across line-intervals par- 
allel, to the X or y-axes, respectively. Within the simple 
model used in this paper, formation of such walls costs 
no energy and therefore they are very stable. The curved 
interfaces seen in Fig. 3a consist of sectors of high-energy 
domain walls across which 0i and "02 change sign simul- 
taneously. The histograms of <I>-values obtained after 
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FIG. 4: Volume fraction, Vord, of an ordered phase > 
1.5) as a function of time after the quench for a series of 
temperatures T/. 



averaging over 10 independent runs for Tj = 0.5465 Tq 
are shown in Fig. 3c for a sequence of time intervals af- 
ter the quench. There, in contrast to Fig. 2c, the initial 
peak at $ = spreads out within the first 100 MCS and 
the maximum at $ ~ 2 starts to grow by 200 MCS to 
reach a considerable height already at 2 • 10^ MCS. The 
$ ~ 2 peak is very narrow, indicating a nearly perfect 
local order in the system. Additional small, sharp peaks 
at multiples of $ = 0.5 may be traced back to the domain 
walls that cross some of the blocs of 4 x 4 x 4 fcc-cells 
used to prepare Fig. 3c. Slow coarsening processes, which 
in this case set in already for t > 2 • lO'^ MCS result in 
considerable sharpening of the main $ ~ 2 peak. 

To better illustrate the differences between the two or- 
dering scenarios described above we show in Fig. 4 how 
Vord: the volume fraction of the ordered phase, grows 
with time for a number of final temperatures. There, we 
counted as ordered all the blocs (of 4 x 4 x 4 fcc-cells each) 
for which <I> > 1.5. The character of Vord{t) changes from 
a slowly increasing function of t for T/ < 0.911 To to a 
sharp, nearly stepwise growth after some incubation time 
for Tf > 0.972 To. 

In the next two sections we analyze in more detail the 
incubation time and the long time coarsening processes. 



tify E{oo) with the energy of a single ordered domain 
of the size of the whole system, equilibrated at Tf. In 
Figure 5a AE{t) per lattice site is shown for a number of 
final temperatures Tf. Initially, the relaxation of AE{t) 
depends on temperature only weakly. Then, after about 
20 MCS, the curves in Fig. 5a split due to a considerable 
slowing down of the relaxation rate with growing temper- 
ature. In turn, for shallow quenches, ordering proceeds 
in three stages, as displayed already in Fig. 2c. First, 
the system remains disordered up to a certain incuba- 
tion time Tine- The metastability of the ordered phase 
is reflected by plateaus in AE{t), which develop near 
Tf ~ Tgp, and extend with increasing Tf. Second, AE(t) 
drops markedly when overcritical nuclei are produced by 
thermal fluctuations. We define Tinc{T) as precisely that 
moment after the quench at which, for a given tempera- 
ture T, AE{t) drops notably (by about 10 percent) below 
its plateau- value. The ordered nuclei formed within the 
incubation time grow against the disordered bulk until, in 
a third stage, ordered regions meet and slow coarsening 
processes set in. 

In addition to AE{t), we also show in Figure 5b and 
5c how the corresponding first moments of the structure 
factors Soi^i^kx") 

kz,t) change with time. To account 
for the anisotropy illustrated in the domain patterns in 
Fig. 3a and b, we consider the two kinds of structure 
factors 



S\\ {k, t) = ^(5i(fc, 0, 0, t) + 52(0, k, 0, t) + 53(0, 0, k, t)), 



(4) 



and 



S±{k,t) = (^1 (0,91, 92, t) 

+S2{qi,0,q2,t) + Sa{qi,q2,0,t)), (5) 

where denotes here the number of {qi , 52) pairs that 
fulfill Qi + q2 =k^, as well as their first moments 



k\\{t) = 



J2k>okS\\{k,t) 
J2k>oS\\{k,t) ' 



(6) 



IV. INCUBATION TIME 

In the preceding section we introduced the incubation 
time Tine as the characteristic time interval between the 
quench and the observed growth of ordered domains in 
the case of metastability, Tf > Tsp. Let us now look more 
precisely at the processes that take place in the system 
within the initial stages of ordering. One relevant quan- 
tity here is the excess energy, AE{t) = E(t) — E{oo), 
defined as the difference between the actual energy of 
the system, E{t), and the energy it would reach after 
complete equilibration, E{oo). This means that we iden- 



k±it) 



'Ek>nkS±{k,t) 

j:k>oS±{k,t) • 



(7) 



In Eq. (6) k denotes one wave-vector component, see 
Eq. (4), while in Eq. (7) k is defined as in Eq. (5). 
Clearly, 5|| {k, t) is determined by structural modulations 
due to high-energy walls, which are reflected, for exam- 
ple, by sign-changes of i/'i when going along the x-axis, 
see Fig. 2. The quantity fc|| therefore characterizes the 
inverse distance between such walls in one symmetry di- 
rection. On the other hand, S± is sensitive to the two- 
dimensional network of low-energy walls, and k± gives 
the inverse distance between such walls, cf. Fig. 3. 
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FIG. 5: Excess energy per lattice site AE (a), and the first 
moments (b) and k±{t) (c) of the structure factors 

Si\{k,t) and S±{k,t) as functions of time after the quench, 
calculated for a series of final temperatures Tf . 



Curves displayed in all three parts of Figure 5 are 
strikingly similar in form, confirming that relaxation pro- 
cesses are indeed dominated by ordered domains growing 
in size. IVIoreover, AE, fc|| and k±, all start to drop be- 
low their plateau-values at the same moment after the 
quench. This fact corroborates our linking the plateaus 
in AE{t) to the temperature-dependent incubation time 
for the formation of overcritical ordered nuclei. In ac- 
cord with the spatial patterns shown in Figs. 2a, 3a and 



3b, where high-energy walls are much further apart from 
each other than the low-energy ones, is in most 

cases considerably smaller than k±{t) calculated at the 
same temperature. 

There are two effects restricting the range of Tinc{Tf) 
accessible to simulations. First, the system size L must 
be considerably larger than the size of overcritical nuclei 
which are formed at temperature Tf. Second, obtain- 
ing smooth curves in Fig. 5 requires a sufficiently large 
number of nucleation events within the maximum com- 
putation time tmax- Again, this becomes increasingly 
difficult to fulfill with growing Tf because both Tine and 
the dispersion of nucleation times over the different sam- 
ples strongly increase with Tf. Comparing our results up 
to tmax = 5 • 10* for lattices of size L=96 and L=128 we 
find Tf ~ 0.973 To to be roughly the highest temperature 
for which Tine can be determined in a reliable way. 

Results of our simulations for Tine{T), plotted in Fig. 6, 
are consistent with the expression 



\nTinc (x{Tq -TfY 



(8) 



which is similar in form to the nucleation rates ob- 
tained within classical nucleation theorjflEa and IVIonte 
Carlo simulations.Eli3 Note that the data point with 
Tf/To = 0.911 falls well below the spinodal, but in Fig. 2 
there still exists a shoulder at that temperature as a rem- 
nant of the plateaus at higher temperatures. 

IncubatioQ. times for CuaAu have been measured by 
Noda et alM from the width F cx k^{t) of the (110) X- 
ray diffraction peak. Close to Tq the width as a function 
of time develops plateaus qualitatively similar to Fig. 5c. 
Their data, however, refer to Tf/T^ > 0.981, whereas 
ours are for Tf/To ^ 0.976. Also one should note that ex- 
perimentally the ratio T^p/Tq may be sample-dependent. 
Because of the sensitivity of Ti„c to temperature differ- 




0.91 0.92 0.93 0.94 0.95 0.96 0.97 0.9 

T.- / rr. 



FIG. 6: Incubation times Tine versus reduced temperature Tf. 
The inset shows the same data (full symbols) after conversion 
to the physical time scale as dessribed in the text, together 
with incubation times measuredEj in CusAu (open symbols). 
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ences with respect to Tgp, this introduces considerable 
uncertainties in our comparison. Nevertheless it is in- 
teresting to present experiments and simulations in one 
plot, as done in the inset of Fig. 6. There the simulation 
data were linked to the physical time scale exactly in the 
manner described towards the end of section II. Both 
sets of data display a remarkable similarity with respect 
to the order of magnitude in time and their trend with 
temperature (and even seem to be a continuation of each 
other). We conclude that the connection between time 
scales for diffusion (see section II) and nucleation as im- 
plied by our algorithm roughly agrees with experiment. 



V. COARSENING REGIME AND SCALING 

As pointed out in the preceding sections, the distances 
between high- and low-energy walls define two differ- 
ent length scales. Moreover, the high-energy walls that 
are much further apart from each other than the low- 
energy ones, are also less stable. As a consequence, we 
expect direction-dependent scaling laws to hold for the 
structure factors Sa{k,t). This is made explicit by con- 
sidering the functions S'n and S± introduced in Eqs. 4 
and 5. Due to their definitions we expect these structure 
factors to obey scaling laws for one-dimensional and two- 
dimensional systems, respectively. This is verified by our 
simulations at Tj — 0.5465 Tq, i. e. in the case of contin- 
uous ordering. As shown by the Id and 2c? scaling plots 
in Figs. 6a and 6b, the data taken for a number of times 
t > 500 MCS after the quench indeed collapse onto single 
master curves. Moreover, the decay of both quantities at 
large k appears to be consistent with Porod's law in d = 1 
and d = 2, respectively ,c 



fc||(05'||(fc,i) (X (fc/fc||(t))- 



ki{t)S^{k,t) (X {k/k^{t)) 



-3 



(9) 



(10) 



We now ask whether the same scaling laws apply to 
the results for Tf > T^p, i. e. to the case of thermally 
activated ordering. Figure 8 contains the corresponding 
plots of S\\{k,t), part (a), and S±{k,t), part (b), cal- 
culated for Tf — 0.972 To at times t > Tine- In both 
parts of the figure scaling regimes are observed which 
are consistent with Eqs. (9) and (10) in an intermediate 
range of fc-values. Both 5|| and S±, however, develop 
tails at large k which show significant deviations from 
scaling. These deviations can be traced back to irregular 
shapes of overcritical nuclei growing against the disor- 
dered bulk (cf. Fig. 2). Such highly structured grain 
boundaries eventually meet, leading to likewise struc- 
tured domain walls when the regime of slow coarsening is 
reached. Small inclusions of disordered phase still remain 
trapped between boundaries of ordered grains for quite 
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FIG. 7: Scaled structure factors k\\S\\{k,t) (a) and k'^S±(k,t) 
(b) calculated for Tf = 0.5465 To < Tsp-, the slopes of the 
dashed straight lines represent a decay according to (fc/fcy)"^ 
(a), and {k/k±)~^ (b). 



a long time after the onset of coarsening processes in the 
bulk of the crystal. 

Closer inspection of our data near the onset of those 
tails reveals that indeed, with increasing time t, the va- 
lidity of Porod's law extends to larger k in both parts of 
Figure 8. We may take this as indication that Porod's 
law, defined by Eqs. (9) and (10), will dominate in the 
late stages of continuous as well as thermally activated 
ordering processes. 

The last point we like to address concerns the ob- 
served growth rates of the characteristic domain-size, 
l{t). Since the transition considered here is described by 
a set of non-conserved order parameters ■01 through ip^, 
one might expect that ordered domains grow accordiag 
to the Lifshitz-AUen-Cahn law l{t) cx t" with j/ = 1/2M 
Experiments on CusAu indeed wpre i.ntfirpreted in terms 
of this conventional growth law.E30'LJ In this context, 
however, one should be aware of the fact that systems of 
the type considered may show effective growth exponents 
during relaxation different from the Lifshitz-Allen-Cahn 
law. In the case of a vacancy mechanism effective expo- 
nents ly > 1/2 in principle can arise in models based-|(pH| 
a non-conserved order parameter. In these studies|2il'E3 
an increased v was ascribed to an accumulation of va- 
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FIG. 8: Same as in Fig. 7 but for Tf = 0.972 To > Tsp. 



cancies within the domain boundaries, leading to an en- 
hancement of the interfacial dynamics. For the present 
model we have verified, however, that the interaction pa- 
rameters chosen do not fayor segregation of vacancies in 
the domain boundaries.Ej Another modification of the 
Lifshitz-AUen-Cahn law that can act in the opposite di- 
rection results from local changes in composition within 
antiphase boundaries so that the ordering process gets 
coupled to a (slow) conserved order parameter compo- 
nent. 

In our simulations we systematically observe effective 
growth expoaents v ^ 1/4 or even smaller (cf. Fig. 5 and 
Kessler et aln). Moreover, within our accessible comput- 
ing times (3 • 10^ — 5 • 10^ MCS), these exponents depend 
slightly on temperature. One possible source of these dif- 
ferences between simulations and the experimental data 
(i/ ~ 1/2) seems to be the existence of the low-energy 
type-I walls. Specifically, let us recall that our model 
contains only nearest neighbor interactions, so that type- 
I walls have zero energy and thus are extremely stable. 
It has been suggestec l pre viously that such a situation 
may lead to v — 1/4.E3E3 Moreover, the present model 
does imply a coupling of the non-conserved order param- 
eters to the conserved density ^pQ . This coupling becomes 
active within type-II walls. The relaxation of a modi- 
fied composition within type-II walls is slowed down fur- 



ther because type-II walls are interconnected via type-I 
walls which have fixed V'o ^'Hd do not allow any exchange 
of composition. The slight upward bending of the low- 
temperature data in Fig. 5 at the longest times, indicat- 
ing an even slower growth, might be interpreted in this 
way, although this point needs to be clarified in further 
studies. 



VI. SUMMARY AND CONCLUSIONS 



Implementing the vacancy mechanism in a model for 
the atomic dynamics in A^B-type fcc-alloys, we investi- 
gated the growth of ordered domains after a temperature 
quench below the transition temperature Tq. Depending 
on the depth of the quench we observe two clearly dis- 
tinct ordering-scenarios: thermally activated nucleation 
of the ordered phase for shallow quenches, Tq > Tf > T^p, 
and spinodal ordering, when Tf < Tgp. Here, the spin- 
odal temperature Tgp was taken over from independent 
simulations at equilibrium. 

In the case of thermally activated processes there is 
some characteristic incubation time. Tine, after which 
a small fraction of the system is covered by overcriti- 
cal nuclei. Detailed simulation results for Tmc, based 
on vacancy-atom exchange, were presented. The time 
period t < Tine manifests itself in plateau regions for en- 
ergy relaxation and for the size of ordered domains, when 
plotted versus t. Clearly, Tine is expected to diverge as 
Tf approaches Tq from below. Within the simple model 
investigated here and the available maximum comput- 
ing time. Tine grows more than 260 times in a narrow 
temperature interval above Tgp - from about 90 MCS at 
Tf = 0.911 To to roughly 24 • lO^MIJS at Tf = 0.976 Tq. 
In comparison with measurementsO of Tine this seems 
to constitute the correct order of magnitude when Monte 
Carlo times are converted to physical times with the help 
of experimental tracer diffusion coefficients. In the coars- 
ening regime we observed growth of the characteristic 
domain-size that was clearly slower than the conven- 
tional t^^^-law expected for curvature driven processes 
in the presence of non-conserved order parameters. This 
may be due to the enhanced stability of the low-energy 
domain walls in the case of nearest-neighbor effective 
interactions assumed in our model. The enhanced stabil- 
ity of low energy walls leads to strong anisotropics of the 
domain-shapes observed in our system, and furthermore, 
to independent scaling laws for the correlation functions 
S'||(fc,t) and S±{k,t). These two functions scale accord- 
ing to Porod's laws for dimensions one and two, respec- 
tively, when Tf ^ Tgp. For Tf > Tsp strong deviations 
from these scaling laws are observed in the region of large 
/c- values, due to highly structured domain walls which 
originate from the stage of fast growth of ordered nuclei 
against disordered bulk. 
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